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Abstract 

QCD coherence effects in initial state radiation at small x in deep inelastic scatter- 
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1 Introduction 



HERA is the first accelerator which allows one to study experimentally the region of low 
values of the Bjorken variable x in deep inelastic lepton-proton scattering {x ~ 10~^ and 
> 1 GeV^). It is, therefore, particularly important to analyze experimental results in 
this kinematical region with the help of theoretical models which contain essential features 
of QCD in the small-x limit. 

Angular ordering (coherence) in initial state gluon radiation (ISR), imposed on real 
and virtual gluon emission at small x, is one of the most important features of pertur- 
bative QCD in the small-x limit [||, |2|, ^. The dominant contribution to ISR is given by 
multi-gluon amplitudes with angular ordering of subsequent gluon emissions. Outside the 
angular ordered regions of phase space destructive interference takes place and amplitudes 
cancel. This phenomenon can be formulated in a probabilistic framework as a branching 
process, in which the gluon structure function at small x is obtained after a summation 
of large log(l/2;) powers corresponding to infrared and collinear singularities coming from 
the angular ordered regions of phase space. The resummation is effectively done with the 
help of the CCFM evolution equation g |]. 

The log(l/x) terms resulting from angular ordering are subleading with respect to 
large logarithms obtained in the multi-Regge kinematics^ leading to the BFKL evolution 
equation Q . In a fully inclusive quantity like the gluon structure function the subleading 
logarithms cancel, and the CCFM and BFKL evolution equation are equivalent. This is not 
true, however, for exclusive processes for which the "angu/ar" logarithms give important 
contribution to the final state structure Q. Therefore, the branching process with angular 
ordering is indispensable for the final state description. 

SMALLX is the Monte Carlo model |5|, |6| which incorporates the CCFM branching 
scheme in the initial state gluon radiation. For an alternative Monte Carlo model which 
also uses the CCFM branching see [^. The dominant process at small x - the quark- 
antiquark pair production accompanied by gluon radiation, see Fig.l, - is generated using 
two approximations for gluon emission. The first approximation, called all-loop and given 
by the CCFM branching scheme with angular ordering for all x values, is compared with a 
conventional branching, called one-loop, leading to the DGLAP evolution equations with 
angular ordering at large x and transverse momenta ordering for small-x values ||8|, ^. 
The two branching schemes coincide at large x but differ significantly when x is small. 
The cross section for the discussed process is obtained using the /c_|_-factorization theorem 
|10| |, in which the off-shell boson-gluon fusion cross section is convoluted with the gluon 
structure function computed in the two branching schemes. 

The aim of this paper is the analysis of angular ordering effects in DIS in the small-x 
HERA kinematics with the help of the SMALLX Monte Carlo program. The comparison 
between the CCFM (all-loop) and DGLAP (one-loop) branching schemes plays an impor- 
tant role in this analysis. The main idea is to fix the parameters of generation in both 
approximations seperately such that a reasonable agreement with the measured inclusive 
cross section is obtained. We choose the charm structure function data from HERA for 

^In this kinematics emitted gluons are strongly ordered in rapidity and have comparable transverse 
momenta. 
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Figure 1: The process generated in SMALLX. 

this purpose. Then, we study various exclusive final state characteristics to show the ef- 
fect of angular ordering. This is done in accord with the present belief that only exclusive 
measurements can discriminate between mechanisms of QCD evolution. 

The paper is organized as follows. In section 2 we briefly describe the process under 
consideration. This includes discussion of the process kinematics, gluon branching schemes 
and their relation to the CCFM and DGLAP evolution equations. In section 3 we describe 
the details of Monte Carlo generation pointing out the crucial elements discussed later in 
section 4 devoted to the results of our Monte Carlo studies. 

2 Description of the process 
2.1 Cross section and kinematics 

The process generated in SMALLX, dominating electron-proton DIS at small value of the 
Bjorken variable xb, is shown in Fig. 1. A quark-antiquark pair is produced accompanied 
by initial state gluon radiation. 

In the high energy limit ^/s ^ M, where M is a quark mass, the cross section for this 
process is given by a convolution of the hard subprocess 'jg — > QQ cross section a, and 
the unintegrated gluon structure function describing gluonic emission |l^] 

a = ^ f ^^^d'k^ {! + (!- yf} .F(x,k^,/.) a{M;k,q) , (1) 
2tt J X y q'^ 
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where y is the standard DIS kinematical variable, q — yl + q± and k = kn — XnP + k±_ are 
the virtual photon and the last exchanged gluon momenta respectively (/ and P are the 
incident lepton and proton momenta), and /U is a hard scale to be discussed later. It is 
important to stress that a is the ofF-shell cross section, which means that the virtuality 
of the gluon, k"^ = k\ = — kj^, is different from zero. This allows a smooth limit of a when 
a quark mass M ^ since the gluon virtulaty k\ regulates a potential divergence present 
in the on-shell version of this cross section, see 1^, ^] for more details. 

The whole information about the initial state gluon emission is contained in the un- 
integrated gluon distribution function k_L, /i) obeying the CCFM equation Q which 
sums angular ordered gluonic emissions, for both large and small xb- This equation allows 
a representation of J- in terms of branching processes, which is convenient for Monte Carlo 
simulation. Before discussing this representation it is necessary to describe kinematics of 
the gluonic emission. 

The momenta of exchanged space- like gluons fcj, and final state real gluons qi = 
ki-i — ki, are decomposed into a longitudinal and transverse components (Sudakov de- 
composition) 

ki = Xip + Xip + k±i , (2) 

where the p and p are null vectors {2p ■ p = s) and k±i ■ p = k±i ■ p = 0. In the frame in 
which p and p are collinear, p ~ (1, 0, 0, 1) and p ~ (1, 0, 0, —1), we have k± = (0, k_L, 0). 

Each emission is characterized by the branching variable Zi < 1, being the fraction of 
the k~^{= ko + ks) component of the exchanged momentum ki-i lost due to emission of 
the gluon qi 

k+ 



Notice the obvious relation 



(4) 



which implies that the small-x limit of the branching process is not necessary related to 
small values of the branching variables. In order to characterize the phase space of multi- 
gluon emission it is convenient to define the rescaled transverse momentum of emitted 
gluons 

q'±i = -rz~ = ^ ^^"^ ' 

where from now q±i = |q_Li| is the length of the transverse component and 9i is the emission 
angle of the massless final state gluon, defined with respect to the z-axis of the collinear 
frame. 

In the approximation leading to the CCFM equation (all-loop approximation) the 
dominant contribution to the multi-gluon emission process comes from a part of phase 
space with angular ordering of subsequent emissions 

Oi > Oi-i ^ q'^i > Zi_i (6) 

where the equivalence results from relations (^). We additionally introduce the collinear 
cutoff go = ■2o9j_o ^^^^ gluon emission: q'^^ > q^. Outside of the angular 
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ordered region destructive interference takes place and multi-gluon emission amplitudes 
cancel. Notice that in the limit Zi <^ 1 condition dH) practically does not restrict the 
ordinary transverse momenta: q±i ~ q'j_- > 0. 

This is not the case in the conventional DGLAP (one-loop) approximation, in which 
a more restrictive condition is imposed 

Qxi > I'u-i : (7) 

and q'^^ > qo. Notice that since q'j_i_i > ^j-i Q'±i-i, condition implies angular ordering 
for all values of z. For ~ 1 conditions (0) and @ are equivalent, but for z ^ the 
transverse momentum ordering q±i > q±i-i is additionally imposed in the one-loop case 
on the angular ordered multi-gluon emission. Therefore, the one-loop phase space is much 
smaller than in the all-loop case. 



2.2 Branching probabilities 

Having the momentum ki-i of the exchanged gluon, the ki and qi momenta can be con- 
structed provided the branching variable and the rescaled tranverse momentum (^) 
are known. Therefore, a chain of subsequent gluon emissions is obtained in which the 
branching variables {zi,q'^-) are generated with the probablity dUi, given in the all-loop 
approximation by 

dUi = As(i, i - 1) Pg{i) Q{q'u - z^-lq'u-l) 0(1 -Zi- Qo/q'u) dz, (8) 

where Ag' and Pg are functions of branching variables in a shorthand notation to be 
discussed below. The first step function imposes angular ordering condition (^) on 
generated real gluon momenta, whereas the second one allows only emissions with the 
transverse momentum q_\_i = (1 — Zi)q'j_^ > Qq to regularize soft singularity in the real 
gluon emission, manifest at Zj = 1 in the function Pg (eq.(p^)). 

The function Ag, called the Sudakov form factor, sums virtual gluon emissions which 
are also angular ordered. This is reflected in the lower limit of the first integration below 

A / / / ^ f f"''' dq'l f^-Qo/q'^ dz } 

^s{q±i, Zi^iq±i^i) = < - ^— / OiSi ■ (9) 

I J(2,_igY^_J2 qf Jo l-z J 

The strong coupling constant as = Sas/ir depends on q± = {1 — z)q'j_, and the upper 
limit of the z integration contains the infrared cutoff Qo which now regularize the soft 
singularity in the virtual emission. The Sudakov form factor has a simple interpretation, 
being the probability of having no gluon radiation within the angular region 9i-i < < 9i. 

The function Pg is the gluon splitting function 

Pgizi,q±i, k^i) = H ANs[Zi, q^^, k^^) , (10) 
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where k±i = \q±i + ... + q±i\ is the total transverse momentum of emitted gluons. The 
function A^s, cahed the non-Sudakov form factor, sums virtual corrections relevant only 
in the small-z limit, with infrared singularities due to soft (zj = 0) exchanged gluons 

^Ns{zi,q\i,k'\_i) = exp ^-asikli) J — J ^0{kj_i - k) 9{k - zq±i) j 

= ..p{-M*L)ln(i)ln(|j-)}. (U, 

where k±i > q±i is assumed to get the last equality. The non-Sudakov form factor screens 
the 1/z singularity in (|lO|) , suppressing radiation both for small z and q±i <C k±i / ^fzl. 

In the one-loop approximation the phase space ordering condition (0) implies that 
Zi-\ = 1 both in the first step function in (^) and in the Sudakov form factor (^). In 
addition, the non-Sudakov form factor A^vs = 1 since the small-z virtual corrections are 
relevant only in the all-loop approximation. 



2.3 Gluon structure function 



The unintegrated gluon structure function k_|_, /i), which obeys the CCFM equation 
1^, § , describes the probablity of finding a gluon with a longitudinal momentum fraction 
X and transverse momentum k^ at the hard scale /U, being related to the maximal angle 
of the emitted gluon. The additional dependence of T on the hard scale is a particular 
feature of the CCFM equation, in contrast to the BFKL equation that, being conformally 
invariant, does not introduce an additional scale into the unintegrated gluon structure 
function. 

The CCFM equation allows one to represent .7-" as a branching process with the help 
of the branching probabilities (§) 

T{x, k±,ij,, qo) = 6{x - xq) 6^{k± - k^o) As{n, qo) ©(/^ - qo) + (12) 



E / / ©(^-^n^ln) ^s(/^,2n9l„) jn^^i^ (5(x-X„)(5^(k_L -kj 
n=l 1=1 



where (x„,k_L„) are parameters of the last exchanged gluon momentum in the chain: 
kn = XnP + kx„, see Fig. 1, and we explicitly indicate the presence of the collinear cutoff 
go for the first real emission. 



The integration in eq. (12) is performed over angular ordered phase space of the real 
gluon emission - the integration variables and their range are implicit in the branching 
probabilities dUi (|8|). The additional G functions introduce a hard scale /i giving the 
upper limit for the phase space integration to be discussed in detail in section 3. The soft 
singularities regularized by the cuttof Qo in dUi cancel after the phase space integration 
and the limit Qo ^ can be taken for !F. The collinear divergence related to the first 
gluon emission and regularized by qo stays as an important parameter in the branching 
procedure, providing the scale for the non-perturbative distribution J'^{xo,i^±o, qo), to be 
convoluted with J^. 
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The structure function T integrated over transverse momentum up to the scale n = Q 
becomes the usual gluon structure function (gluon density) 

Fix,Q,qo) = J d2kx^(x,kx,g,go) e(Q-|kx|) . (13) 

We are interested in the dependence (evolution) of F at small x given by the CCFM 
equation, where for simplicity the fixed strong coupling constant as is assumed. The result, 
given in terms of the moments F^^ of the gluon structure function F, is the following 



FUQ,Qo) = dx x^ F{x,Q,qo) = i^^j , (14) 



where the gluon anomalous dimension 7 is equal to that obtained in the analysis of the 
BFKL equation 



,(^)=7— (^) = ^ + 2C3(^)' + 2C5(^)% (15) 

and Ci is the Riemann zeta function. 

It should be noted that the BFKL anomalous dimension was derived in the all-loop 
(CCFM) approximation as a result of the complete cancellation of large log(l/x) powers 
coming from the integration over the angular ordered phase space of the real and virtual 
emissions. These logarithms give subleading corrections to the BFKL equation and are 
fully cancelled in the gluon structure function at small x. They are not cancelled, however, 
for more exclusive quantities, and in such cases the all-loop approximation with angular 
ordering is essential for the proper description of exclusive processes, see Q for recent 
discussion. The subleading corrections in the CCFM equation also decrease the BFKL 
value of the QCD Pomeron intercept loq = (4 In 2)0:5 in the asymptotic formula: J-'{x) ~ 
^-(i+i^o)^ and reduce diffusion in the transverse momentum space, see [11, O. 



In summary, the CCFM equation is a generalization of the BFKL equation by tak- 
ing into account a part of subleading log(l/x) corrections generated by angular ordering. 
These corrections cancel in the fully inclusive quantity like the gluon structure function 
( [T3| ) but give an important contribution to more exclusive quantities (e.g. gluon multi- 
plicity). 

In the one-loop approximation the small-x behaviour is given by the gluon anomalous 
dimension equal to the first term in expansion (p!5[): 7(0;^/^) = a^/w- This result is 
in a perfect agreement with the gluon anomalous dimension obtained from the DGLAP 
equation in the small- a; limit. 



3 Monte Carlo simulation 

The cross section (|l]) and branching formulae ( [T2|) form the basis for the Monte Carlo 
generation in SMALLX. 

The procedure starts with generation of an initial state parton shower, see Fig.l. In 
the first step the gluon momentum /cq = xoP + kj_Q is generated according to a given input 
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^(xo,k_Lo). Then next gluon emissions are generated with distributions (^). At the i*'' 
vertex, is selected with the help of the Sudakov form factor @, and an azimuthal 
angle is chosen randomly to form a vector q^j- The branching fraction Zi is chosen next 
according to distribution (^). Notice that Pg depends on the variable k^i = |k^j_i + q±j| 
thus, in order to simplify the generation, the approximation q^j ~ q^j is made when the 
selection of Zi is done. Then the emitted gluon momentum q^j = (1 — Zi)c(^^ is computed 
and the true value of k_|_j calculated. 

The branching continues until the stopping condition is reached, i.e. for certain n we 
have 

Znq'±n = Xn Vs tan(6'„/2) < /i , (16) 

and q'±n+i momentum in the next emission exeeds the bound. In this case the n-gluon 
final state is generated. In the all-loop approximation the bound fi results from condition 
® 

9i < ... <e„ <e , (17) 

where is the emission angle of the quark-antiquark total momentum pq + pQ = Yp + 
Yp + Q±. Thus, it is easy to see that the choice 

^ = XnVs tan(e/2) (18) 



gives the angular ordering condition (17). In the one- loop approximation the stopping 
condition takes the form 



= V^' + M2 , (19) 

where ^ is given by the scale at which the hard cross section a starts to be strongly 
suppressed as a function of the gluon virtuality k'j_. We will show in the forthcoming 
analysis that condition (^) is much more restrictive than (16) and leads to smaller values 
of transverse momenta than in the all-loop case. 

The quark-antiquark pair can be produced provided 

s = {k + qf > , (20) 

where k = kn is momentum of the last gluon in the chain and q yl + q± is the photon 
momentum, see Fig. 1, generated with the distribution dy/y dq\/q\. If condition ( pOD is 
satisfied the hard cross section a in eq. (|l|) is calculated and the event is accepted with a 
non-zero weight, otherwise the event weight is equal to zero. 

In general, the non-zero weight is a product of subweights coming from the virtual 
photon state generation, the gluon branching and the hard scattering. At each stage 
of generation the subweights could become zero when momenta fall outside kinematical 
limits. Condition (|20| ) is quite restrictive in this respect and only a fraction of gluonic 
events can produce a quark-antiquark pair. 

At the end of the generation the longitudinal momentum carried by gluons is computed 
at the scale /x 

P(^) = I dx xF{x,p.) , (21) 







where F is given by eq. ([T^). The conservation of P(/u) is slightly violated when p, runs, 
thus eventually all weights are renormalized by the factor 0.5/ P{p) to get gluons carrying 
half of the proton's momentum. 
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4 Monte Carlo results 



We shall present the results of our first experience with the SMALLX Monte Carlo program 
in HERA kinematics (27.6 GeV + 820 GeV for electron and proton energy, respectively). 
Most of the results will be presented in the HERA laboratory frame in which the z-axis 
is defined by the incoming proton direction^. Some results are presented in the virtual 
photon-proton frame (CMS frame) with the z-axis also defined by the proton direction. 

We confine ourselves to the study of open charm production in DIS at small x, with the 
charm quark mass M = 1.5 GeV. This is motivated by our strategy to fix the generation 
parameters, in both approximations separately, in order to describe the measured at HERA 
inclusive quantities like structure functions. Then we compare more exclusive final state 
characteristics to understand the difference between the two gluon branching schemes. 
Within such an approach open charm production is the best process to consider since it 
is not spoiled by the details of the light fiavours generation. 

We fix the generation parameters to describe the charm contribution to the pro- 
ton structure function F2 measured in DIS processes at HERA This includes the 
parameters of the initial gluon distribution , for which the following form is postulated 

xo.FO(xo, k^o, Qs) = ^0 V (1 - xoY' exp(-2kio/Q?) , (22) 

where the parameter Qg also plays the role of the collinear cutoff go for the first real gluon 
emission, and the normalization constant Aq is determined from the momentum sum rule 
P{Qs) = 0.5 (|2l|). We have fixed = 1 GeV and pi = 5 motivated by large-x behaviour 
of the conventional gluon distribution. For the purpose of Monte Carlo generation we 
keep the non-zero value of the infrared parameter Qq in ^ giving the minimal transverse 
momenta of emitted gluons: Qq = 1 GeV. The parameter pQ is chosen to obtain a good 
description of at = 12 GeV^, whereas the result for higher values is given by 
the evolution embodied in both approximations. In order to achieve that, the shapes of 
the initial distributions have to differ significantly: for the one-loop mode the small-x 
values are strongly enhanced (po = 0.4) in comparison to the all-loop initial distribution 
{po = 0.0), see the corresponding lower curves in Fig|2[ 

In Fig. I we show F^^^ at = 12, 25 and 45 GeV^ for all-loop (soUd fines) and 
one-loop (dotted lines) approximations with the parameters defined above. The evolution 
does not differ significantly once the parameters have been tuned at = 12 GeV^, and 
both modes describe ^2^^ in the whole measured range quite well. This result indicates 
that there is no significant difference between the evolution in the two approximations 
when the asymptotic x-behaviour of is obtained. Notice that in the one-loop approx- 
imation the asymptotic behaviour is imposed by hand, through the steep initial gluon 
distribution, whereas in the all-loop mode is a result of the CCFM evolution - compare 
the two corresponding upper curves in Fig. |2[ The dashed lines in Fig. |3|, obtained from 
the flat [pq = 0.0) input distribution, illustrate the importance of the steep form of the 
initial distribution in the one-loop mode for the data description. 

^ This frame is also used in SMALLX to define decomposition (jij): p = P and p — I, where P and I are 
the colhding proton and lepton momenta, respectively. 
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4.1 Gluon branching properties 



Having fixed the generation parameters we compare different characteristics of the gluon 
branching in the two approximations for two cases. In the first case, we analyze the gluonic 
state generated until the stopping condition (^) or ( p^ ) is satisfied, with the event weights 
given only by the gluon branching procedure. In the second case, we take into account 
only those gluonic events which satifisfy condition (|2^) for quark pair production, and 
the full event weight, with the 75 — > QQ subprocess contribution cj, is computed. We 
analysed the small-x DIS regime, 1.6 • 10-"^ <XB < 1.3 • 10-2 for < Q2 >= 12 GeV^ and 
0.01 < y < 0.8, based on the statistics of 10^ generated events. 

In Fig. 1^ we show the distribution of the transverse and longitudinal components of the 
last exchanged gluon momentum /c„ = x„P + k±, see Fig. 1. The most striking feature of 
transverse momentum distribution (two upper plots) is a dramatic difference in the hard 
part of the spectrum (k± > 5 GeV), where the hard tail in the all-loop approximation is 
strongly enhanced (solid lines); see also analysis ([^). This phenomenon can be traced 
back to the all-loop condition ([T7| ) which imposes only an angular bound on the maximal 
angle for emitted gluons, in contrast to the one- loop condition ( p!9| ) that also bounds 
their transverse momenta. Since kj_ is the total transverse momentum of emitted gluons, 
k± = \q±i + ■•. + Q±n\, this effect is also clearly visible in the presented spectra. There is 
no difference in the longitudinal momentum spectra between the two approximation (two 
lower plots) as a result of the choice of the input distributions. 

The events with the charm pair produced ("with bo3^' plots) have suppressed small 
values of Xn and large values of k± in comparison to the "no 6ox" plots, where all generated 
gluonic events are analysed. This can be easily understood if condition (20) is written for 
q ^ yl + q± and k x„P + k± 

Clearly, large values of Xn and small values of |k^| are favoured by quark-antiquark pair 
production. Therefore, the gluon distribution is dynamically suppressed in the generation 
of the full final state for x = Xn < 10^^, as is seen in the lower "wii/i 602^' plot in Fig. 
see also 

The analysis of the final state real gluons, done in the same spirit, is presented in 
Fig. |. The two upper plots with the gluon multiplicity distribution exhibit a pronounced 
difference between the branching schemes. Due to angular ordering (|6|) there are signif- 
icantly more gluons in the all-loop scheme, and their transverse momenta can be much 
bigger than in the one-loop case - see the plots in the second row, where the hard tail in 
the all-loop g_L-spectrum persists. Notice that most gluons are "soft" {q± ~ 1 GeV) and 
there are significantly more "soft" gluons in the all-loop approximation. The sharp edge 
at 1 GeV is the result of the infrared cutoff Qo imposed on transverse momenta. The plots 
in the bottom row show transverse momentum as a function of rapidity {E± flow). The 
difference between the two approximations in these plots reflects the fact that there are 
more gluons (both "soft" and "hard") in the CCFM branching. 

It is interesting to analyse the role of the hard scale (18) in the all- loop approximation, 
related to the maximal angle for the gluon emission. In Fig. I we show the same charac- 
teristics as in Fig. I, but now the one-loop stopping condition ( [T9| ) is used in the all-loop 
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generation. The g_L-spectrum reveals that the scale ( [l8| ) is responsible for the hard tail 
in the all-loop approximation - the hard tail disappears when the one-loop bound (|l^ is 
imposed in the all-loop mode, while the "soft" part of the spectrum stays intact. The other 
all-loop characteristics are not so dramatically affected. The average number of gluons 
is smaller and there are less radiation events with more than n = 4 gluons in the "wii/i 
6oaf' case (upper plots). The significant difference between the generation modes persists 
in the E± flow (bottom plots) since the "soft" part of the g_L-spectrum is not affected by 
the one-loop stopping condition. The lack of the "hard" gluons {q± > 5GeV) influences 
mostly the E± flow in the central rapidity region — 1 < ?] < 1. In summary, angular or- 
dering ( p!?! ) in the all-loop approximation leads to significantly more gluons radiated into 
the final state with the A;j_-spectrum much harder than in the one-loop approximation. 



4.2 Final state analysis 

In this section we study the full final state properties, including the charm-anticharm 
quark pair. Therefore, for each analysed event condition (^) is satisfied and the final 
weight is computed. Before presenting the results let us stress that we study the very 
nature of the CCFM and DGLAP branching schemes, since the additional effects like the 
final state radiation from quarks and gluons or hadronization effects are neglected. These 
effects remain to be included in SMALLX to make the comparison with the real exclusive 
data. 

In Fig. we show two final state characteristics discussed in the context of the BFKL 



signatures in the HERA data: the transverse momentum spectra [0, 15 1 and the tran- 
vserse energy flow [^], for the small-x kinematics specified in the previous section. The 
presented plots reveal how important the charm and anticharm quarks are for the final 
state characteristics. The difference between the approximations in the hard part of the 
gluonic g_L-spectrum is significantly reduced by the quark g^-distribution. However, the 
difference between the approximations persists, favouring harder spectrum in the all-loop 
mode which stays in a qualitative agreement with the recent HERA data ||l^. The E± 
flow is less affected by the quark inclusion, especially in the forward region, but we expect 
signiflcant hadronization corrections in this region. 

In Fig. 1^ we show the mean transverse energy flow (Ej-) in the central region of 
pseudorapidity in the 7p CMS frame for several values of x. (Ej-) decreases with increasing 
X for the all-loop generation mode, and increases in the one-loop case. This result is in a 



qualitative agreement with the semi-analytical calculations done in [16|, which are based 
on the BFKL and DGLAP equations. 

The diffusion of the gluon transverse momenta k± along the gluon chain as Xn decreases 
can result in large values of k_i_ in the all-loop approximation, see Fig. ^. This phenomenon 
leads to azimuthal decorrelation between the quark and antiquark momenta . In Fig. ^ 
we show the distribution of the difference A(j) between the azimuthal angles of the quark 
and antiquark momenta, calculated in the 7p CMS frame, for two different values of 
the Bjorken variable xb- In the all-loop mode the distribution shows weakening of the 
azimuthal (back-to-back) correlation as x decreases, while for the one- loop case this effect 
is not observed, and the distribution is concentrated around = tt, independent of x 
values. 
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In Fig. 10 we present another widely advertised signature of the BFKL dynamics at 
HERA: the cross section for forward jet production [^], where the jets were chosen as in 
the HI analysis [|^. In this process a forward jet with respect to the proton direction 
with transverse momentum g_L ~ Q is measured, where Q is the virtuality of the photon. 
The cross section for such a measurement is significantly enhanced as a result of gluonic 
emissions populating the rapidity region between the forward jet and the virtual photon. 
The original proposition |jl^ concerns the gluonic emission in the Regge kinematics, but 
in view of discussion in the angular ordered kinematics imposed in the CCFM scheme 
is more appropriate. We remind the reader that all our results refer to charm production, 
so that the presented forward jet cross sections are much smaller than those presented in 

in 



The total forward jet cross sections, shown in the upper plot of Fig. 10 as a function of 
xb for the two generation modes, differ significantly. The all-loop cross section is bigger 
and much steeper than in the one-loop approximation. The two lower plots with the pure 
gluonic and pure quark forward jet cross sections explain the reason for the difference seen 
in the upper plot. In the all-loop approximation the total cross section is dominated by the 
gluonic forward jet contribution and the pure quark component is almost negligible for the 
smallest xb values, whereas in the one- loop mode the gluon and quark contributions are 
comparable, and much smaller than the gluonic cross section in the all-loop case. Clearly, 
forward jet production is a good process to discriminate between the two gluon branching 
schemes. 



5 Conclusions and outlook 

We have presented Monte Carlo studies of angular ordering effects in initial state gluon 
radiation in the small-x HERA kinematics, described by the CCFM (all-loop) branching 
scheme. The studies have shown the relevance of the CCFM approximation in the de- 
scription of exclusive processes in DIS at HERA. In contrast to the conventional DGLAP 
(one-loop) approximation the CCFM scheme allows more gluons in the final state with 
a much broader transverse momentum spectrum. As we have shown this is the result of 
the angular ordering condition ( p!7|) imposed on the phase space of gluon emissions in the 
all-loop approximation. 

The gluonic spectra are altered by the quark pair production but the angular ordering 
effects are still visible in the presented characteristics. The most spectacular results are the 
azimuthal angle decorrelation for the heavy quarks - a direct consequence of A;_|_-diffusion 
along the gluon chain in the CCFM scheme, and forward jet production. The rest of the 
results may strongly depend on such effects like the final state branching, hadronization or 
inclusion of quarks in initial state radiation. Therefore the additional effort is necessary to 
include these effects into SMALLX, before the meaningful comparison with HERA data 
could be done. We hope to report on progress in this direction in the near future. 
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Figure 2: The integrated gluon distribution F in all- loop (solid) and one-loop (dotted) 
schemes at the initial scale 1 GeV^ (the corresponding lower curves), and after evolution 
to < >= 12 GeV^ (the corresponding upper curves). 
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Figure 3: ^2^^ in the all-loop (solid lines) and one-loop (dotted lines) approximations. The 
dashed lines are obtained for the flat input distribution in the one-loop mode. 
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Figure 4: The transverse k± and longitudinal Xn momentum distributions of the last 
exchanged gluon for all gluonic events ("no 6oaf' ) and the events with the quark pair 
produced {"with boo^^) in the all-loop (solid) and one-loop (dotted) approximations. 
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Figure 5: The final state gluon characteristics in the all- loop (solid) and one-loop (dotted) 
approximations. The influence of the quark pair production is shown on the right ("wii/i 
6oaf' plots). 
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Figure 6: The same as Fig. ^ but now the one-loop stopping condition is imposed in the 
all-loop generation. 
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Figure 7: The A;_L-spectrum and Et flow for final state partons. Gluon and quark contri- 
butions are shown separately. 
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Figure 8: The average value of the transverse energy in the central rapidity region —0.5 < 
r]* < 0.5 in the jp CMS frame as a function of the Bjorken variable x. 
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Figure 9: The difference between the azimuthal angles of the quark and antiquark momenta 
in the jp CMS frame for different values of the Bjorken variable x. 
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Figure 10: The forward jet cross section as a function of the Bjorken variable xb, with the 
HI analysis cuts. The gluon and quark forward jet cross sections are shown on separate 
plots. 
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